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Abstract. Consider a semiclassical Hamiltonian 

H v>h := h 2 A + V - E 

where h > is a semiclassical parameter, A is the positive Laplacian on K d , V £ C'£°(M. d ), 
i.e. V is a smooth, compactly supported potential function which is central, that is, 
depends only on \x\ and E > is an energy level. In this setting the scattering matrix 
Sh(E) can be defined at any positive energy E; it is a unitary operator on L 2 (S d— 1 ), 
hence with spectrum lying on the unit circle; moreover, the spectrum is discrete except 
at 1. 

We show under certain additional assumptions on the potential that the eigenvalues 
of Sh(E) can be divided into two classes: a finite number ~ c ( j(R\/E/h) d ~ 1 , as h — > 0, 
where B(0, R) is the convex hull of the support of the potential, that equidistribute 
around the unit circle, and the remainder that are all very close to 1. Semiclassically, 
these are related to the rays that meet the support of, and hence are scattered by, the 
potential, and those that do not meet the support of the potential, respectively. 

A similar property is shown for the obstacle problem in the case that the obstacle is 
the ball of radius R. 



1. Introduction 

In this paper we consider the scattering matrix for a semiclassical potential scattering 
problem with spherical symmetry on K d , d > 2. Let V be a smooth, compactly supported 
potential function which is central, i.e. V(x) depends only on We consider the Hamil- 
tonian 

(1.1) H Vth := h 2 A + V - E 

where A = — Yli=i ®f * s ^ e positive Laplacian on R d , E > is a positive constant (energy) 
and h > is a semiclassical parameter. 

The scattering matrix Sh(E) for this Hamiltonian is defined in terms of the asymptotics 
of generalized eigenfunctions of Hv,h- For each function qi n 6 C 00 (S C ' _1 ), there is a unique 
solution to Hy ^u = of the form 

(1.2) u = r-^' 2 (e-^' h qM + e + ^' h q out {-uj)) + 0(r^ d+ ^ 2 ), 

as r -> oo, see e.g. [15] . Here q out € C 00 ^ -1 ). The map q in i-> e'^^^^oat is by defini- 
tion the scattering matrix Sh(E). The factor e 17 ^- 1 )/ 2 j s chosen so that this 'stationary' 
definition agrees with time-dependent definitions (see e.g. [T7] or [Hj), and is such that the 
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scattering matrix for the potential V = is the identity. It is standard that the scattering 
matrix Sh(E) is, for every h > 0, a unitary operator on L 2 (§ rf_1 ), and that Sh{E) — Id is 
compact. It follows that the spectrum lies on the unit circle, consists only of eigenvalues, 
and is discrete except at 1. It is therefore possible to count the number of eigenvalues of 
Sh(E) in any closed interval of the unit circle not containing 1. In fact, semiclassically (i.e. 
as h — > 0) we are able to separate the spectrum of Sh(E) into two parts. One is associ- 
ated to the rays that meet the support of the potential; to leading order in h there are 
Cd{RVE J 'h)^ 1 of these, Cd = 2/{d— 1)!, and the other part is associated to the rays that 
do not meet the support of the potential. Those eigenvalues corresponding to rays that do 
not meet the support are close to 1, as one should expect, since the eigenvalues of the zero 



potential are all 1 — see Proposition 1.3 The other eigenvalues are affected by the potential, 
and we can ask whether these 'nontrivial' eigenvalues are asymptotically equidistributed on 
the unit circle. Indeed Steve Zelditch posed this question to one of the authors several years 
ago. 

Before stating the main result, we discuss further the scattering matrix in the case of 
central potentials. In this case the eigenfunctions of the scattering matrix are spherical 
harmonics and the Helmholtz equation reduces to an ODE in the radial variable. The 
generalized eigenfunctions then take the form u = r~( d_2 " 2 /(r)17™, where 



(i.3) -a 2 --^- ^,^ +iv±-z\ f = 




and /(r) = Hj^, d _ 2)/2 (ry/E/h) + c(l)H^, d _ 2)/2 (rVE/h) for r > R, where the i# are the 
standard Hankel functions, [T]. With our normalization, ShY™ = c(/)y" ; m . In particular, the 
eigenvalue o/K m is independent ofm. We write the eigenvalue corresponding to Y" ; m in the 
form e ll3l ' h . The quantities Pi,h/2 are frequently referred to as 'phase shifts.' See e.g. [T7] 
for a review of these facts. 

We now discuss the conditions we on the potentials in the main theorem. Both conditions 
are dynamical conditions, i.e. conditions on the Hamiltonian dynamical system determined 
by the symbol of Hyh- As usual in microlocal analysis we refer to the classical trajectories 
of this system as bicharacteristics. We first define the interaction region 

(1.4) TZ := {x : V(\y\) < E for all \y\ > \x\} . 

This is the region of x-space accessible by bicharacteristics coming from infinity. Notice that 
for central potentials this region takes the form 

(1.5) TZ = {\x\ > r } where r = mi{s > r V(s) < E}. 

r>0 

The first condition is 

(1.6) V is nontrapping at energy E in the interaction region. 

That is, x tends to infinity along every bicharacteristic in TZ both forwards and backwards 
in time. 

The second concerns the scattering angle. Let R be such that B(0, R) is the smallest ball 
containing the support of V, i.e. 

(1.7) -8(0, R) = chsuppV, the convex hull of the support of V. 

We recall (see Section [2] for further details) that for a central potential, the scattering angle 
E(a) is a function only of the angular momentum a and measures the difference between 
the incident and final directions of the trajectory (which are well-defined, since the motion 
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is free for |x| > R where B(0, R) is the support of the potential) — see [T7]. The scattering 
angle is zero for all trajectories with a > R. Our second condition is that 



(1.8) 

Then our main result is 



The zeroes of S'fa) = — E(a) are finite in [0, R). 
da 



Theorem 1.1. Let R be as in (1.7 1, and assume that V G C™{R d ) be central and sat- 
isfies conditions (1.6) and (1.8 1. Then as h \. 0, we consider the eigenvalues e l ^ l ' h for 
which I < R/h, counted with multiplicity Pd(l) = dimker (Agi-i — 1(1 + d — 2)). There are 
2(Ry r E/h) d - 1 /(d-l)\ + 0(h-( d - 2 ')) of these, and they equidistribute around the unit circle, 
meaning that 



(1.9) 



sup 

O<0 o <0i<2ir 



2(RVE/h) d - 1 /{d-l)l 



2tt 



as h 4- 0, 



where N(<f)o,<f)i) is the number of (3i t h with I < R/h and 4>q < /3i t h < tf>i (mod 2it), counted 
with multiplicity. 



Remark 1.2. Many potentials satisfy conditions (1.6| and |L8j). For example, take any 
smooth compactly supported function W(r) such that W(r) = for r > R, W(r) > 
for < r < R and W"(r) is positive and monotone decreasing in some nonempty interval 
[R— e,R)- Then cW(r) satisfies the conditions for sufficiently large c. An explicit example 
is 

V/(- 2 -« 2 ), r <R 
0, r>R. 



W{r) 



More generally, in Section [5] we show that a potential satisfying V' < and V + rV" > on 
the interaction region 1Z also satisfies conditions (1.6 1 and (1.8). This condition is actually 



much stronger than ( |1.8[ ): it implies that S'(a) has no zeroes at all on [0,R). 
We also show 

Proposition 1.3. Let k € (0, 1). The eigenvalues e l/3 ' h for I > (R + h K )/h satisfy 



1 



0{h c 



h^O. 



Here and below, 0(h°°) denotes a quantity that is bounded by C^h N for all N and some 
C N > 0. 

Note that the number of eigenvalues not covered by the Theorem and Proposition above 
is o(h 1 ~ d ), and hence cannot affect the equidistribution properties. Hence we get 



Corollary 1.4. Suppose that V satisfies conditions (1.6) and (1. 

we have 



Then for each e > 0, 



(1.10) 



sup 

e<0o<0i^27r— e 



2(RVE/h) d - 1 /(d-l)\ 



2?r 



as h 4- 0, 



where N(<j)o,(/)i) is the number of fti t h (with no condition on I), counted with multiplicity, 
satisfying <pg < (i\^ < <f>\ (mod 2ir). 
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Results directly analogous to those for semiclassical potentials are also true in the case 
of scattering by a disk of radius R centered at the origin. The scattering matrix in this case 
can be defined similarly; given any function gi n £ C°°(S d_1 ), there is a unique solution u to 
the equation (A — fc 2 ) u = such thaiQ 

u = r-^ 1 '/ 2 (e- ikr q in + e+ ikr q out ) + 0(r)-( d+1 V 2 , r -+ oo 

u\\x\=R = 0. 

The scattering matrix Sk is again defined qi n i-> e I7r ^" 1 ^ 2 g out , and the standard facts 
about the operator Sh also hold for Sk- As above, the spherical harmonics diagonalize the 
scattering matrix. We write SkY" 1 = e lXl - k Y [ m . We will prove 

Theorem 1.5. As k —¥ oo, the eigenvalues e lXl - k for which I < Rk, counted with multiplicity 
Pd(l) = dimker (A§d-i — 1(1 + d — 2)), equidistribute around S 1 in the sense of the Main 
Theorem. In fact, we have the stronger statement 

N(4> ,4>i) (l)i-4 



(1.12) sup 

O<</>o<0i<2tt 



2(Rk) d - 1 /{d-l)l 2tt 



0(k~ 1/3 ) as k 



As far as we are aware, the present paper is the first in the mathematical literature to deal 
with the question of the equidistribution of phase shifts over the unit circle. However, there 
are a number of previous studies of high-energy or semiclassical asymptotics of eigenvalues 
of the scattering matrix. The relation between the sojourn time and high-frequency asymp- 
totics of the scattering matrix was observed in classical papers by Guillemin and Majda. 
Melrose and Zworski [TB] showed that for fixed h > the absolute scattering matrix for a 
Schrodinger operator on a scattering, or asymptotically conic, manifold is an FIO associ- 
ated to the geodesic flow on the manifold at infinity for time tt. Alexandrova [2] studied 
the scattering matrix for a nontrapping semiclassical Schrodinger operator, and showed that 
localized to finite frequency, it is a semiclassical FIO associated to the limiting Hamilton 
flow relation at infinity, which includes the behavior of the Hamilton flow in compact sets. A 
more global description was given in Hassell-Wunsch [9] where the semiclassical asymptotics 
of the scattering matrix were unified with the singularities of the scattering matrix at fixed 
frequency (i.e. the Melrose-Zworski result [16]). These results are explained in Sections [2] 
and [3] below. 

Asymptotics of phase shifts, i.e. the eigenvalues of the scattering matrix, were analysed 
by Birman-Yafaev [21 [H [51 IB] , Yafaev [TH] and more recently Bulger-Pushnitski [7 . This 
was done in the context of the scattering theory for a fixed potential. In this case, the 
scattering matrix S(X) tends in operator norm to the identity as A — s- oo so the phase shifts 
tend to zero uniformly. The asymptotics of the phase shifts for a fixed energy, and also the 
high-energy asymptotics, were analyzed. 

In |21) . Zelditch and Zworski analyzed the pair correlation function for eigenvalues of 
the scattering matrix associated to a rotationally invariant surface with a conic singularity 
and cylindrical ends. They showed that a full measure set of a 2-parameter family of such 
surfaces obeyed Poisson statistics, as conjectured by Smilansky, motivated by the Berry- 
Tabor conjecture for completely integrable Hamiltonian systems. 

In a slightly different setting Zelditch |20j analyzed quantized contact transformations, 
which are families of unitary maps on finite dimensional spaces with dimension N — > oo. 
He proved under the assumption that the set of periodic points of the transformation has 



^Here we prefer to use non-semiclassical notation where the energy level is fc 2 , as is traditional in obstacle 
scattering literature. The variable fc here corresponds to 1/h above, when the energy level E = 1. 
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measure zero, that the eigenvalues of these unitary operators becomes equidistributed as 
N — > oo. After reading a draft of the current paper, Zelditch pointed out to the authors 
that a similar strategy could be used in the context of semiclassical potential scattering to 
prove equidistribution. In fact, this strategy is likely to be a more direct approach to proving 
equidistribution than the one we employ here. On the other hand, our approach has several 
advantages: it also gives approximations to the individual phase shifts, up to an 0(h) error 



see Proposition 3.1 ), and in addition appears to be a better method for obtaining a rate of 



equidistribution, as in Theorem 1.5 above 



2. Dynamics 

We now review some standard material on Hamiltonian dynamics for central potentials. 
Consider first the case the dimension d = 2. 

The classical Hamiltonian corresponding to our quantum system is 

\e + V(r)-E 

or in polar coordinates, using (r, (p) and dual coordinates (p 7 77) , 

2 

H = p 2 + \ + V{r) - E. 

The Hamilton equations of motion are 

r = 2p 

(2-1) r 2 

f) = 0. 

The invariance of the Hamiltonian under rotations is reflected in the conservation of angular 
momentum r\ — 2r 2 (p. For a given bicharacteristic, this is the minimum value of r along the 
free bicharacteristic that agrees with the given one as t — > — oo (we could just as well take 
t — » +00 since it is a conserved quantity). 

Notice that in the case of general dimension d, each bicharacteristic lies entirely in a 
two-dimensional subspace, so the above discussion in fact includes the general case. 

The scattering matrix is related to the asymptotic properties of the bicharacteristic flow. 
Geometrically this information is contained in a submanifold 

L C T'S* -1 x TPS* -1 x R 

that we define now. Returning to the case of general dimension d, we identify § d_1 with the 
unit sphere in M. d and identify the cotangent space T* § d_1 with the orthogonal hyperplane 
uj 1 - to lu. Given lu and rj £ T^S d_1 , take the unique bicharacteristic ray whose projection 
Xu>,ri(t) to IR'' is given by rj + tui for t << 0. Define (uj' , rf) by 

uj (oj, rj) — lim x(t)/\x(t)\ 

t— loo 

rj (uj,rf) = lim x(t) — (x(t),uj')uj' 



(2-2) 



and t(uj, rj) to be the sojourn time or time delay for 7; this is by definition the limit 
(2.3) lim ti(a) — t^a) — 2a = t(uj, rj), 

a— >oo,a>>0 



(i 
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Figure 1. Here x u>v is the classical trajectory equal to rj + tuj for t « 0. 
The scattering angle £(77) is the difference between the outgoing direction 
uj' and the incoming direction uj. Note that \rj\ = \rj'\ by conservation of 
angular momentum since the potential V is central, V = V(r). 



where ti(a) is the smallest time, t, for which r(t) = a and £2(0) is the largest. We then 
define L to be the submanifold 

(2.4) L := {(u, n, u/(u, rj), -rf(u, T}),r).} 

As shown in [9], L is a Legendrian submanifold of T* § d_1 x T* § d_1 x K with respect to 
the contact form x + x' — dr, where \ is the standard contact form on T* S d_1 , given in 
any local coordinates x and dual coordinates £ by £ • cfa;. Note that the projection of L 
to T* S^" 1 x T* S d_1 is Lagrangian with respect to the standard symplectic form. Indeed 
it is the graph of a symplectic transformation {u),rf) <— > (ui'^r)'), and the scattering matrix 
is a semiclassical Fourier Integral Operator associated to this symplectic graph [2], [S]. 
The sojourn time, however, carries extra information and is directly related to high-energy 
scattering asymptotics as observed in [14], [8], [9]. 
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The previous paragraph applies to any potential, central or not. We return to the case of a 
central potential, for which, as observed above, the dynamics take place in a two-dimensional 
subspace, so we can assume d — 2 without loss of generality. In that case we use the angular 
variables <p, ip' in dimension d = 2 instead of w, ui' above. Consider a bicharacteristic with 
angular momentum 77 £ R and initial direction ip e S 1 . The scattering angle £(77) is 
defined to be the angle between the initial and final directions, normalized so that S is 
continuous and £(77) = for 77 > i?, i.e. 

(2.5) £(77) = <p'(<p, 77) - ip, £(77) = for 77 > R, 
(see Figure [T]) and the sojourn time is defined by 

(2.6) T( V )=r(p, V ). 



Notice that both £ and T depend only on 77 in the central case. The fact that L in (2.4) is 
Legcndrian then implies the following relation between these functions: 

(2.7) dr = r, ■ (dtp - dtp') =► -^T(ry) = -t^Efo). 

drj dr] 

Remark 2.1. Notice that the ambiguity of £ modulo 2ir is eliminated by our convention that 
£(77) = for 77 > R. We point out that by reflection symmetry, we have £(77) = —£(—77) 
modulo 2n, but it might not be the case that £(77) = —£(—77) on the nose: this will happen 
if and only if £(0) = 0, which will be the case if and only if the interaction region is the 
whole of Mr. However, we always have £'(77) = £'(—7/), which shows that T'(rf) is an odd 
function, and hence T{rf) is even in 77. 

3. AN ASYMPTOTIC FORMULA FOR THE EIGENVALUES OF Sh 

In this section we find an asymptotic formula for the eigenvalues e 4 ^ ! fe of Sh in terms of 
the scattering angle £. 

Proposition 3.1. Define the real-valued function G(a), a€l, by 

dC 

(3.1) — (a)=£(a), G(a) = for a > R, 

da 



where £ is the scattering angle function in (2.5). Then we have the following estimate on 
each eigenvalue e l ^'- h of Sj 



h ■ 



(i) If the dimension d is even, then for any e > there exists C = C(e, d) such that, for 
all I satisfying e < Ih < R, we have an estimate 



(3.2) 



exp{i( G ((, + ^ W )} 



< Ch. 



(ii) If the dimension d > 2 is odd, then for any e > there exists C = C(e,d) such that 



(3.2) holds whenever a = Ih > e is distance at least e from the set 



(3.3) {a : £(a) € {nk} keZ } . 

To prove the proposition, we will use the fact, proven in [9], [2], that the integral kernel 
of Sh is an oscillatory integral associated (in a manner we describe directly) to the Legendrc 



submanifold L in (2.4). To be precise, the Schwartz kernel of Sh can be decomposed following 



[SJ Prop. 15] (with minor changes in notation) as 

S h =Ki+K 2 + K 3 , 

with the Ki as follows. 
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Fix i?2 > Ri > R- First, K2 is a pseudodifferential operator of order zero (both in the 
sense of semiclassical order and differential order), microsupported in {\ri\ > Ri}, hence 
taking the form in local coordinates z on § d_1 

(2wh)-^-^ [ e l( - z - z ^</ h b{z,(,h)d( 



for some smooth symbol b(z, (, h) equal to zero for |Clg(z) < Ri where | • | g ( 2 ) is the norm on 
T*S rf_1 . This reflects the fact that the Legendrian submanifold L in (2.4) is the diagonal 



relation w = w',77 = —t/,t = for \r)\, \rf\ > R, to which pseudodifferential operators are 
associated. Moreover, K2 is microlocally equal to the identity for \r/\ > R2, i.e. b = l+0(h°°) 
for |C| S > i?2. Indeed, the full symbol (up to 0(/i°°)) of the scattering matrix is determined 
by transport equations along the rays with \rj\ > R. Since these transport equations are 
identical to those for the zero potential, the scattering matrix in this microlocal region is 
microlocally identical to that for the zero potential, which is the identity operator. 

Next, K\ is a semiclassical Fourier integral operator of semiclassical order with compact 
microsupport in {\rj\ < -R2}. That is, K± is given by a sum of terms taking the form in local 
coordinates 

(3.4) Ki(w,u/,h) = ft-^" 1 )/ 2 -^/ 2 f e l ^^'^/ h a(u;,uj',v,h)dv 

Jr n 

with respect to a suitable phase function $ and smooth compactly supported function a. 
Here the phase function parametrizes L locally, meaning 

(1) On the set Crit$ := {(u),u',v) : D v $(u,w',v) = 0}, A^w>$ has rank N. This 
implies that 

(3.5) £(*) := (u,D u ${u,<J,v),(J,D^*(u,rf,v),Q(u,<J,v)) 

is a smooth submanifold. 

(2) L($) n suppa = L n suppa. 

By K\ having compact microsupport in the set {\rj\ < R2}, we mean specifically that if 
(u,r],u)',r/,T) £ L has \r/\ = \rj\' > R 2 and (u,lj',v) € Crit($) with {D u>u >$(u, U)' , vj, $) = 
(77, t]',t), then a{ui, oj', v, h) — 0(h°°) in a neighbourhood of (w, ui',v). 

Finally, K3 is a kernel in C 00 ^^ 1 x § d_1 x [0, ho)), i.e. smooth and vanishing to all 
orders at h = 0. 

For the proof of Proposition |3.1| we need to know the principal symbol of K\ . 



Lemma 3.2. The Maslov bundle of the canonical relation C of the FIO Ki is canonically 
trivial, and with respect to this canonical trivialization, the principal symbol of K\ is equal 
to 1, as a multiple of the Liouville half- density on C coming from either the left or right 
projection of C to T* S d_1 . That is to say, 

(3.6) cr(A'i)(w,?7,w',-?/) = \du:dr]duj' dr]'\ 1/2 , 

for (oj, r], u)', — if) € C such that \rj\ < R%. 

Proof. Consider first the Maslov bundle of C. Notice that C is almost the same as L; in 
fact, it is given by 

C = {(v, r\,J , — ?/) I 3 r such that (w, r), w', rj ', t) e L}. 
Since C is a canonical graph (i.e. the graph of a symplectomorphism), associated to the 



scattering relation as in (2.2 1, it projects diffeomorphically to T* § d 1 via both the left 
and right projections, and the lift of Liouville measure on T* § d_1 via the left projection 
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agrees with the lift via the right projection (since C is a Lagrangian submanifold of T* § d_1 x 
T* and the Liouville measure can be expressed in terms of the symplectic form T* § d_1 ), 
providing a canonical half-density on C. We also note that the scattering relation is the 
identity whenever \r]\ > R since then the corresponding bicharacteristic is not affected by 
the potential. Therefore, over this part of C there is a canonical trivialization of the Maslov 
bundle. Since the Maslov bundle is flat, we can use parallel transport to extend this to a 
global trivialization: in fact, in the case d = 2, the space T* S 1 retracts to S^ 1 x {ij > R}, 
while for d > 3, T* S d_1 is simply connected, hence in either case parallel transport provides 
an unambiguous trivialization. 

We now consider the principal symbol of the scattering matrix. The scattering matrix 
may be viewed as a 'boundary value' (after removing a vanishing factor and an oscillatory 
term) of the Poisson operator, as in [9_, Section 7.7 and Section 15]. The principal symbol 
of the scattering matrix is correspondingly derived from the principal symbol of the Poisson 
operator. The principal symbol of the Poisson operator is real: it solves a real transport 
equation with initial condition 1. Therefore, the principal symbol of the scattering matrix 
is real, up to Maslov factors, i.e. it is a real number times an eighth root of unity. On the 
other hand, unitarity of the scattering matrix shows that the principal symbol lies on the 
unit circle (as a multiple of the canonical half-density); hence it is an eighth root of unity. 
Finally, the principal symbol of the scattering matrix is equal to 1 for \rj\ > R, since here the 
scattering matrix is microlocally equal to the scattering matrix for the zero potential, which 
is certainly equal to 1. Since the principal symbol is smooth, is restricted to eighth roots of 
unity, and is 1 for \r)\ > R, it follows that the principal symbol is equal to 1 everywhere. □ 



Proof of Proposition 3. 1 : First we reduce the problem to the cases d — 2 and d = 3 as 



follows. Writing pi t h,d for the eigenvalue j3i t h in dimension d, observe that by (1.3), 

(3.7) Pi,hA+2k = Pi+k,h,d for d > 2, fc > 0. 

It follows that for d > 4 even, we have fii,h,d = A+(d-2)/2./i,2 and for d > 5 odd, we have 

Pl,h,d = A+(d-3)/2,h,3- 

Consider the case dimension d = 2. For any smooth function G : R — > K, the function 

(3.8) $(<p,<p',v) = ( l p-<p')v + G(v), 



locally parametrizes (see (3.5)) the Legendrian 



dG dG 

:= {(ip,r],(p',r)',T) \ <q = v = -r/,(p' - <p = ^r( v ), T = ~ v ^ v ^ v "> +G ( V )}- 
With G as in JO, this gives an explicit global parametrization of the Legendrian subman- 



ifold L in (2.4) if we take <p € [0, 2n], ip' £ R. In this case the relation between r and E 



given by the last equation above is 



-tj£(»7) + G{ji) => -j-r = -rj— Efa), 



in agreement with (2.7). Therefore, plugging (3.8) into (3.4), the operator K\ takes the 
form 

Ki{<p,tfih) = {2nh)- 1 ( e i( y {[f, ^' )v+G{v) ) /h a{ip-ip' ,v,h)dv, <p G [0, 2n], tp' e K, 



where a is smooth and supported in \v\ < i?2- Notice that we may assume that a depends 
only on (ip — ip',v,h) since the scattering matrix and the phase function both have this 
property. 
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Now we obtain an expression for the eigenvalue e l ^ uh of the scattering matrix Sh on 
Yi = (2ir)~ 1 ' 2 e u ' p using the pairing 

(3.9) e^' h = (S h Y u Yi) = (K^Yfi + (K 2 Y U Y{) + (K Z Y U Y{). 

Clearly (K^Y^Yi) = 0(h°°). Consider the K\ term. Writing I = a/h gives 

(K\Yi,Yi) = (27r/i) _1 (27r) _1 f / / jlMWHi-V))/^ _ v> dv d(p d(p '_ 



o 



Changing integration variables to (ip, (p — tp — ip'), the kernel is independent of the first of 
these variables, so that integrating in it simply removes the factor 2tt. We are left with 

{KiY h Yi) = {2nh)- 1 ( ( e r ^ v+G ^~ a ^' h a{(p,v,h)dvdp. 



The phase is stationary at the point v = a, —(p = G'{v) = S(f) and the stationary phase 
lemma shows that the integral is equal to 

(3.10) {K{Y u Yi) = e tGia)/h a(-Z(a), a, 0) + 0(h) 

(noting that the Hessian of the phase function has determinant 1 and signature 0). 
Next we write 



Here, the phase is stationary when a — v. However, b is supported where |u| > Ri > R while 
a < R by hypothesis, so there are no stationary points on the support of the integrand. It 



follows that (K 2 Y h Yi) = 0(h°°). Thus by (3.9| 



(3.11) e ll3 '- h = e lG(Q)//l a(-E(a), a, 0) + 0(h). 

The principal symbol of K\ as an FIO is given as a multiple of the Liouville half-density 
on T* §\ Idpdrjl 1 '/ 2 , by [H3 Section 3] 

(3.12) a{Kx)(<p, V,<P + Av),~V) = a(~^(v),V, 0)\dpd v \^ 2 . 

Indeed, the density dc defined on page 143 of that paper equals \dipdrj\, where we used 
coordinates (x, 9) — (ip, ip', v). The principal symbol is the image of the map from C to A 
defined immediately following the definition of dc, in the notation of that paper. In the 
notation of the current paper, C = Crit($) and A is the projection of the Legendrian L 



onto the first four coordinates. It follows from equation (3.121 and equation (3.6) that 
(3.13) o(-S(» J ),» 7 ,0) = l. 



Combining (3.13) and (3.10), we see that 

(3.14) e^ l - h =e lG{a)/h + 0(h), 



establishing (3.2) 



We proceed to the case d — 3. In this case, we will obtain the eigenvalue e l ^ l - h by 
pairing the scattering matrix Sh with the highest weight spherical harmonics Y l 1 . These 
concentrate along a great circle 7, which we parametrize by arclength, ip G [0, 2tt]. Choose 
Euclidean coordinates in M 3, so that the two-plane spanned by 7 is the plane x$ = 0. Then 
Y"/ = Ci(xi + ix 2 ) 1 where ci is a normalization factor, equal to (27r)~ 1 / 2 (7rZ) 1 / 4 (l + 0(l -1 )). 
Let 9 be the spherical coordinate equal to the angle with the positive X3 axis. Then we can 
write 

(3.15) Yt(<p, 9) = cie il,p (smey = Cl e Ulp e' l9[e) 
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where g(9) = — log sin 9 



7r/2) 2 /2 + O((0-7r/2) 4 ). 



In particular, expression (3.15) shows (and it is in any case well known) that the Yj' 



concentrate semiclassically (i.e. as I — > oo) at the set {9 — 0, £ = 0, a = a} where / = 
a/h + 0(1). Here we use coordinates (a, dual to (ip., 9). So, to compute the pairing (3.9) 



with Yf replacing Yj, we first need to determine an oscillatory integral expression for K\ 
that is valid in this microfocal region. (Note that the K% and K3 terms give an 0(h°°) 
contribution as before.) So choose «o distance > e from the set (3.3). As we will see, it 
suffices to find a local parametrization of L in a neighbourhood of 

{9 = 9' = 7r/2, ( = C' = 0, <t = -a' = a , ip - ip 1 = S(a )}; 

this is the set of incoming and outgoing data of bicharacteristics which remain in the x% = 
plane. To define this parametrization, we consider first a parametrization in two dimensions 
locally near a bicharacteristic where r\ — a. As we have seen such a two dimensional 
parametrization is (ip— ip')v+G(v), for v close to a. We note that when v = a, ip' — ip = E(a), 
and we can write it in the form 



(3.16) 



ip — ip = ± dist(<y9, ip ) + 2irk 



for some integer k (recalling that the distance dist((y9, ip') lies strictly between and ir). We 
now claim that a suitable phase function is 



(3.17) 



$(uj,lj',v) = (=F dist(w, u') - 2nk)v + G(v), 



where v G R is localized near ao, G(v) is as in (3.1), and the sign =p and the value of 
k agree with the two-dimensional case. Indeed, on each two-plane, if we use spherical 
coordinates ('ip, 9) adapted to that 2-plane then the form of the phase function agrees by 
construction with the two-dimensional phase function and therefore parametrizes that part 
of L associated to that 2-plane (since the dynamics on each 2-plane is identical to the d = 2 
dynamics), that is, the subset (in the coordinates adapted to that 2-plane, indicated by a 
bar) 



(3.18) 



{9 = 9 =7r/2, C = C =0,^ -Xp = Y,(a),a 



T(a)}. 



We now observe that we can eliminate k by redefining G(v) locally to be G(v) + 2nkv, 
which only has the irrelevant effect of changing E by 2irk (notice also that th is d oes not 
affect the eigenvalue formula involving e lG ( lh )/ h m the statement of Proposition 3.1 ). From 
here on we only work with the 
that this means that < ip' — 



sign in (3.16), i.e. the 



sign in (|3.17|), and k = 0. Notice 
Returning to 



ip < 7T and < E(a) < 7r, i.e. sinE(a) > 0. 
our spherical coordinates associated to the 2-plane X3 = 0, we can use the spherical cosine 
law applied to the spherical triangle with vertices (ip, 9), (ip 1 , 9'), and the pole X3 = 1: 

cos dist((cp, 9), (<p' , 9')) = (cos(ip — tp') sin 9 sin#' + cos 9 cos 9') 



to write 

(3.19) $((/3, y, tp , y, v) — — cos -1 (cos(<,£> — tp') sin 9 sin 9' + cos 9 cos 9') v + G(v). 
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We can then write in these coordinates 



L = 



a = = — — — - 

sin cast (w, u) ) 



( = dg<$> 



-r (sin(ip — ip') sin sin fl 7 ) 



(3.20) 



sindist(u>,a/) 



— (cos(p> — ip') cos 9 sin 6' — sin 9 cos 9') 



a' = cU$ 



— — - — : "7 fsinfw — tp') sin 9 sin 6') 

sindist(w,u/) v VY ^ r ' ' 



C = dg.Q = - 



— (cos(tp — tp') cos 9' sin 9 — sin 9' cos ( 



sin dist(w, u' 

t = dist(tj, ui')v + G(v)\ where dist(w,a;') = G'(v). 

(Notice that by direct inspection we see that this agrees with (3.18) when 9 = 9'= tt/2, 
since then cos 9 = cos 9' = and dist(w, u/) = ip' — ip = S(a) and so a = v = a.) 

The scattering matrix, microlocalized to this region of phase space, will then take the 
form 



(3.21) 



(2irhy 3 / 2 / e i *( w ' w ''">/ h o(w,w',»,/i)di;. 



In terms of this parametrization the principal symbol of (3.21 ), say where both u> and io' 

d(ip,9,a,(,d v $) -i/2 



lie near the great circle 7 and hence where we can use coordinates (ip, 8, ip' ' , 9 '; a, £, c', C', t), 
is given at the point (ip, tt/2, ip + S(a)/i, 7r/2, a, 0, —a, 0, r(a)) by [TU] 



(3.22) a(cp, n/2, tp + S(a), tt/2, a, 0)e~" /4 Ids d<9 da dC 



1/2 



det 



where the e~ l7r / 4 is a Maslov factor; see Remark 3.3 for more discussion about this. We 
need to compute the determinant above. We can disregard the repeated coordinates (<p, 9) 
and compute, using (3.20), 
(3.23) 

0-1 




det a if' C ^? =det| 



sin(< y 3 / — <p) 

1 G"(v)j 



d(tp',9',v) 
It follows that the principal symbol is 

(3.24) a(ip, tt/2, ip + S(a), tt/2, a, 0)e~ i7r/4 ( 
Then by equation ( |3.6| 

(3.25) a(<p, tt/2, 93 + S(a), tt/2, a, 0)e" i7r/4 = 



sm(ip' — (p) sinE(o;) 



at 9 = 9' = 



vsin S(a) 



-V2, 



c?s d# dcr dQ 



1/2 



f — 1 
VsinEfaW 



1/2 



sin E(a) - 

We next write the contribution of K\ to the expression (13. 91) for the eigenvalue e^' ! 



Writing I = a/h and using (3.15) we get 
(3.26) 

(iri^y/) - m- 3 / 2 y e i*M,<p',o , ,v)/h e -ia( V - V ')/h ^y /2 (27r )-i 
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Here the factors (a/7r/i) 1 / 2 (27r) _1 are to normalize the functions Yf in L 2 . We will analyze 
this using the stationary phase lemma with complex phase function, see e.g. (Ill thm. 7.7.5]. 
Here the phase is 



(3.27) 



9, tp', 9', v)=$-a(tp- tp') + ia(g(9) + g(9 r )). 



Notice that the integral as a function of (tp, tp') depends only on <p — tp' by the rotational 
invariance of the scattering matrix, and the form of the Y/ which take the form e il(p times 
a function of 9. We change variable to (tp,tp), tp = ip — tp' and integrate out the variable 
tp, giving us a factor of 2tt. Then \E' has nondegenerate stationary points in the remaining 
variables (tp,9 ,9' ,v). The imaginary part of the phase is stationary only at 6 = 6' = tt/2, 
while stationarity of the real part requires that v — a and — tp = G'(v) = T,(a). The 
stationary phase lemma then gives us that (3.26) is equal to 

\2 



(3.28) 



2tt (Wr 3/2 (^) 1/2 c^r 1 ) ( 2 ^) 2 



1 



det(-iD 2 4') 1 /: 



: a(tp, tt/2, tp + E(a), vr/2, a, 0) + Q{h) 



Here, to keep track of constants, we have written out all constants in (3.26); the first 2tt 
comes from the integral in tp and the (2ttK) 2 comes from the leading term in stationary 
phase in the four variables (tp,9,6' ,v). Simplifying the constants and using (3.25) this is 
equal to 

V2 / ^G(a)/h \ 

0(h)). 



(3.29) 



2a 2 



sinE(a) 



(det -iD 2 ^) 1 / 2 



We will show that, in the above expression 
(3.30) det-iD 2 ^(ip,0,tp + Y l (a),0,a) 



2ia 2 
sin Ti(a) 



iE(c 



Accepting this for the moment, we obtain from (3.28) 

e m, h = e iGla)/h e iE(a)/2 + Q^y 

Since = G'(a) and a = Ih, this can be written 



(3.31) 



JGUl+l/2)h)/h 



0(h) 



completing the proof of Proposition 3.1 



It remains to prove the formula for the Hessian in (3.30). First we notice that when 
9 = 9' = tt/2 we have, using the formula in ( 3.19 ) for the distance function, in the coordinates 

&,e,e'), 



o 



o 



(3.32) 



L> 2 dist(( / 3, 7 r/2,^ + E(a),7r/2) = cotE(a) 





esc E(a) 



— cscE(a) cotE (a) 



From this (3.17) and (3.27), we conclude that in the (v,tp,y,y') coordinates 



(3.33) 



D ^ = 



/ G"(v) 1 

1 

— acotT,(a) + ia 

\ acscE(a) 






a esc E(a) 
-a cot E(a) + ia 
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Thus 



det -iD 2 ^ = -a 2 ( cot 2 E(a) - 2i cot S(a) 



E(a) 



2za 2 



sin £(a) 



-iS(o) 



and (3.30) holds. 



□ 



Remark 3.3. The Maslov factor in (3.22) and (3.25) arises as follows. First, Lemma 3.2 



shows that the Maslov bundle over L is canonically trivial. However, unlike in the case 
d = 2, there is a nontrivial Maslov factor from comparing our phase function $ above to one 
— let us call it $ — that agrees with the canonical phase function, i.e. the pseudodifferential 
phase function, for \r)\ > R. By [TDJ Theorem 3.2.1], the principal symbol written relative 
to $ contains the Maslov factor e l7r<T / 4 where a is the difference of signatures, 

a = sgnL> 2 t ,$-sgnD^$ 

where w = (wi,W2) are the phase variables for A tedious computation shows that 
a = — 1, leading to the Maslov factor in (3.22) and (3.25). (We remark that since $ depends 
on one phase variable and $ on two phase variables, by [101 Equation (3.2.12)] a is odd, so 
the Maslov factor cannot vanish in this case.) Of course, the Maslov factors are irrelevant 
to the question of equidistribution, but they are relevant to the question of determining the 
eigenvalues modulo 0(h). 

Proof of Proposition \1.3[ In view of the remarks in the proof of Proposition |3.1[ specifically 
equation (3.7), it is only necessary to do this in the cases d = 2 and d = 3. For dehniteness, 
we write down the proof for d = 3; it is similar, and in fact simpler, for d — 2. Consider 
a spherical harmonic F/ with hi > R + h K , where n < 1. The eigenvalue e % ^ l - h is given by 
pi?] ) with F/ replacing Y\ . 

First assume that hi > R' > R. Then the K\ term in (3.9) will be 0(h°°) (for a suitable 
decomposition of Sh = K\ + + K$ as above, with R^ < R'), so we only have to consider 
the K2 term. This is given by a pseudodifferential operator with symbol equal to l + 0(h°°), 
so the (Y l , K2Y1 ) term is equal to 1 + 0(h°°), proving the Proposition in this case. 

Next assume that R + h K < hi < R[. For R' < R u the K 2 term in (|3j]) will be 0(h°°) 
(for some other decomposition of Sh, with Ri > R'), so we only need to consider the K\ 
term. That is, it remains to show that 

((Ki - Id)F/, y/) = 0(h°°) for R + h K < hi 

Using as above polar coordinates (f,6) on S 2 with dual coordinate (a, Q, we find a phase 
function ^ for K\ that parametrizes L microlocally in the region 

{\(a,0\ g >R-S} 

for fixed small S > 0. Indeed, since L is given by the diagonal relation 

(3.34) { (p = p>,e = 6',o- = -a',C = -C,T = 0} for {\(a, Q\ g > R} , 

it follows that the functions (<p,8,<j' , £') furnish local coordinates on the Legendrian L for 
C) Is > R} an d therefore, by continuity, for {|(er, C)| g > R — 8} for some small S > 0. It 
then follows from [T2l Theorem 21.2.18] that L can be parametrized by a phase function of 
the form 

— ip v — 9 w + H(ip, 9, v, w). 

Since K x is pseudodifferential for \(a' , (')\ 2 g = C' 2 + (sin 9')- 2 a' 2 > R 2 (i.e. L satisfies pl4| |), 
we have 

v>Rsm9' =^ H = ipv + 9w and b = 1 + 0(h°°) 
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Thus 
(3.35) 



(2tt/i)3 + tA j ' 



As above we have written I — a/h; hence a > R + h K . 
We insert cutoff functions by writing 

fv-Rsm6'\ , .fv-RsmO' 
1= X rz. +(1-X) 



h K ) \ hK 

where \(t) is supported in t < 1/2, equal to 1 for t < 1/4. With the cutoff x inserted, the 
phase function is nonstationary on the support of the integrand, since stationarity requires 
that v = a. It follows that we can integrate by parts arbitrarily many times, using the fact 
that the differential operator 

1 h d 

v — a i dip' 

leaves both exponential factors invariant; doing this gains a factor of h 1 ~ K each time since 
a — v > h K /2 on the support of the integrand. Thus the \ term is 0(h°°). 
With the cutoff 1 — x inserted, we write the integral in (3.351 in the form 



(3.36) 



j(^( v -ip')v+(e-0')w)/h^i(^H(^,e,v,w)-ipv-ew)/hf ) ^/ qi v w ^ _ 

x C J\ v ~ Rs[n9 ' \ p i a {g{e)+g{e'))/h dipdedy'de' dvdw 

l{ X '\ h K J (2tt/i)3 



We claim that the factor 



j( H { v ,e,v,w)- lp v-ew)/h b ^ (p ^ ^ Wj Wj ^ 



1 x(l-x) 



v - R sin 6 1 



is 0(h°°). In fact, the term in the large brackets is 0(h°°) for v > i?sin(9', while if v < 
RsinO', then the 1 — x term vanishes identically. It follows that the 1 — x term is also 
0(h°°), completing the proof of Proposition 1.3 □ 



4. EQUIDISTRIBUTION 



For clarity of exposition, we will assume in this section that the energy E in ( 1.1 1 satisfies 

E = l. 

The lemmas and propositions in this section apply to the arbitrary E > case by writing 
V = V/E and h = h/VE, 



If u = {e 
defined by 

(4.1) 



. ,e 



2-jtixk 



} is any set of K points on S 1 , then the discrepancy D(w) is 



D(lj) :— sup 

0<(t> <rj> 1 <2Tr 



N((j>o,(j)i;u)) 0i-0o 



K 



2tt 



where N((f>Q,4>i;uj) is the number of points in uj with argument in [0q,0i] (modulo 2tt), 
counted with multiplicity. We state the following lemma in slightly more generality then is 
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necessary for semiclassical potentials so that we may apply it without significant modification 
to the case of scattering by the disk. 

Lemma 4.1. Let G(a) : [0, R) — > R be smooth and assume that 

{a : G (a) — 0} is finite in [0, R). 
Consider the points xik on the unit circle 

(4.2) E h = {x lk := exp (i(G(lh)/h) : < Ih < R, k = 1, . . . ,p d (l)} , 

included according to multiplicity. Here Pd(l) — dimker (A§d-i — 1(1 + d — 2)). 
Then the sets Eh equidistribute as h J, 0. That is, the discrepancy satisfies 

(4.3) lim D(S h ) = 0. 

To apply the lemma to the eigenvalues of the scattering matrix Sh , we must show that they 
still equidistribute despite satisfying only the weaker asymptotic condition in Proposition 
|3.1| Below, S will denote any set satisfying 

(1) S is a finite subset of [0,R]. 

(2) If < a < R and G"(a) = 0, then ae S. 

(3) 0,ReS. 



For example, to apply Proposition 4.2 to Proposition 3.1 when the dimension d is odd one 
takes 5* to be the set to be 

{a < R : G"{a) = 0} U {a < R : G'{a) = 0} U {0, R} . 

Proposition 4.2. Notation and assumptions as in Lemma \4.1\ let 

(4.4) £ h = {x lk :0<lh<R,k = l,... , Pd (l)} , 

be a collection of points on S 1 (included according to multiplicity), such that for any e > 0, 
if I satisfies dist(Z/i, S) > e then 

xi k = expi(G(lh)/h) + E(l, h)) 



where \E(l,h)\ < C(e)h. Then, if G satisfies the condition in Lemma 1^.1 



Xm\D(E yi ) = 0. 

We will use the following notation. With any set S as above, let 

E h (e) := £ h n {x lk : dist(//i, S) > e} 

(4.5) 

£ h (e) :=E h n{xi k : dist{lh,S) > e}, 
always understood to include points according to multiplicity. 



Proof of Proposition ^. 6 ^ assuming Lemma \4-.1\ By Proposition [3TTJ for every e > 0, there is 
a constant C = C(e, S) > so that for h sufficiently small, 

N{fo + Ch, 0i - Ch\ Sh(e)) < N(<t> , 4>r,£~h{e)) < N{fo - Ch, fa + Ch; £ h (e)). 
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More precisely, C depends on the number of points in S H [e, R — e]. Dividing through by 
2(R/h) d - 1 /(d - 1)! and subtracting (fa - fa)/2n gives 



N(<t> Q ,fa;£ h (e)) fa-<po 



2(R/h) d - 1 /(d~l)\ 



2tt 



< max 



N(<p + Ch,fa-Ch;£ h (e)) fa - fa 



2{R/h) d - 1 /(d-l)l 
N(cf> -Ch,fa+Ch;£ h (e)) 



2/r 



2(R/h) d - 1 /(d-l)l 
< D(£ h (e)) + (1 + C(e))0(h) + 0(e), 

uniformly in h and e, where for the second inequality we used 

2R d ~ 1 



2rr 



(4.6) 



|4(e)l = ^(ri-l)! (1 + Q(ft) + Q(e)) ' 



where |£zi(e)| is the number of points in £/t(e). Similarly D(£h) — D(£h(e)) + 0(h) + O(e) 
for h, e small. The same is true for £/,. Thus 



(4.7) 
Thus 



N(fa,fa;£ h ) 
2(R/h) d ~ 1 /(d-l)\ 



>i - <Po 
2ir 



<D(£ h ) + (l + C(e))0(h)+0(e). 



limsup_D(£ft) = 0(e), 



and as e > was arbitrary, we obtain the result. 



□ 



Remark 4.3. Note that the proof gives no information about the exact vanishing rate of 
D(£h) as h —¥ 0. For this, one must have information on the dependence of C(e) on e, and 
then optimize in e in (4.7 1 as h — > 0. This is what we do in Section [6] to obtain improved 
remainders in the case of scattering by the disk. 

To prove the lemma, we use theorems from [13] . The following theorem follows from [13, 
ch. 2, eq. 2.42] 

Theorem 4.4 (Erdos-Turan). There is a constant c > such that if 

uj = {e 27rlXl ,...,e 27riXN } 
is a finite sequence of N points on S 1 and m is any positive integer, then 

(4.8) D(u) 



To bound the exponential sums that appear on the right hand side of (4 
ch. 1, thm. 2.7], namely 



Theorem 4.5. Let a and b be integers with a < b, and let f be twice differentiable on [a, b] 
with \f"(x)\ > p > for x £ [a,b]. Then 



(4.9) 



if{l) 



<(|/'(6)-/'(o)|+2)(— +3 



Finally there is [El thm. 2.6] (with minor modifications in notation) 
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Theorem 4.6. For 1 < i < k, let uji be a set of points on S 1 with discrepancy D{iOi). 
Let uj be a concatenation of U)\, . . . , that is, a set obtained by listing in some order the 
terms of the CJ{ . Then 

k 

(4.10) D(w) < £ HD( Ui ), 

i=l ' ' 

where \ui\ is the number of points in u>. 



Proof of Lemma J^.l We begin by assuming that G" has no zeroes in the open interval 
(0, R). In this case, the set S in (4.5) above is equal to {0,R}. 
We first analyze the subset £/j(e) C £h defined in Define 

P = P(e) 

(4.H) „ \ 



mm \G"[a)\ 

e<a<R—e 

2 max \G'{a) 

e<a<R—e 



We will show that for each 7 G (0, 1) there is a constant c = 0(7) > so that for each e > 0, 

(4.12) D(8 h (e)) <c(jP + Kp-^h 1 ' 2 ^' 2 + p- 1/2 h 1 ' 2 + hh 1 ^) 

Since £fr(e) — = /i~ d+1 (0(e) + 0(h)), for some c = 0(7) > independent of e we have. 



(4.13) 
So 



D(£ h ) 



sup 

O<0 o <0i<2tt 



N((f>o,(f>i;£h) 4>\-4>g 



\£h\ 

<c(e + h)+D(£ h (e)) 



2tt 



lim sup £>(£/! ) < ce, 



and as e > was arbitrary, we obtain (4.3). 

Case 1: dimension d = 2. Note that when d = 2 the multiplicity of the eigenspaces is 
P2(l) = 1 if ' = and 2 otherwise, so that 

|f h (e)|=2(L(J2-e)/ftJ-re//»l+l). 



We apply Theorem |4.4| with uj — £/j(e), so that, in the notation of Theorem 4.4 
G(lh)/{2irh). Thus 



x, 



rn 

D(£ h (e))<c[-+Y / ^ 



1 



l(R-e)/h] 



1 E 



JjG(lh)/h 



_(«-,)AJ-V/*l + . ,„ 

Then we apply Theorem [45] with /(at) = (j/2ir)G(xh)/h, a = \e/h~\, and &=[_(#- e)//ij. 
Thus, if x/i < i? — e then = /17 |G"(ir/i)| /27r > hjp/2n, which equals p in the notation 

of Theorem 14.51 It follows that 



(4.14) 



D{£ h (e)) < c 




By letting m = \_h 7 J for any 7 > 0, we obtain (4.12). 



Finally, suppose there are a finite number of points < 01 < . . . < o,n-i < R with 
G(ai = 0, and let ao = 0, a n = R. Note that, if we define £h(a,b) to be the set of 
Xi t k with a < Ih < b, counted with multiplicity. Then the above arguments show that 
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lim/j D(£h{a, b)) = 0, in fact if p(e, h) (resp. k) is defined to be min Q+e < Q ,<b_ e |G"(a)| (resp. 
min), then the proof is the same. The lemma in the d = 2 case now follows from Theorem 4.6 
since by this equation 

n 

D{£ h ) <Y / D ( £ h(a l -i,a i )). 



The proof is now complete in the case d = 2. 

Case 2: dimension d > 2. As in the d = 2 case, we begin by assuming that G"{a) has 
zeroes only at and R. 



We will apply Theorem 4.6 to D{£h{e)) decomposed as a superposition in the following 
way. It will be convenient to set 



(4.15) 
Define, 



N:= [(R-e)/h\. 



with unit multiplicity. Note that u)(n) has N — n + 1 elements. Setting 

oji = w(0) 
uj 2 = w(l) 

uj pd(1) =u(l) 
w P<i(i)+i = w ( 2 ) 

Up d (2) = w(2) 
w Pd(W-l) + l = v( N ) 



w Pd(A) = w(iV), 

we see that the set £/i(e) is the superposition of the sets wi, . . . u> Pd (N)- 

The discrepency D(u>(n)) can be estimated using the method from the d = 2 case. In 
particular, as in (4.14) we see that for any positive integer m, 



(4.16) 



D(ui(n)) < c 



y * 



1/2 
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By Theorem 4.6 we have 



Pd(N) 



i=0 



A 



(4.17) 



^ E i^^i -p<*( n - 



n=0 
N 



< ch d - x ^(N -n+ l)(n + l) d - 3 L>(w(n)), 



71=1 



Substituting the estimate (4.16) into (4.171, again with m = \h~^\ for some fixed 7 <G (0, 1), 
we end up with five terms to deal with corresponding to the five terms in the right hand side 
of (4.16). For all of these we use standard bounds for sums of polynomials and N ~ c/h. 
The easiest is the 1/m term, since 



N 



(4.18) 



h d ~ l J2(N - n + l)(n + l) d ~ 3 /i 7 < c/1 7 . 



Next we do the terms involving p. There is 



(4.19) 



and 



A 



[h->] 



/i d " 1 ^(7V-n + l)(n+l) d - 3 



J=l 



(N-n + l)j {Kj) {jcph) 



( 32 \ 



1/2 



n=l j=l WP 7 

/ 1 x 1/2 A \h—>\ 

n=l j=l 



1/2 



1/2 



< c 



ph / 

/1 \ 1/2 

(i) »-^/», 



A 



(4.20) 



The other terms are 



h d - 1 Y > (N -n + l)(n + l) d " 3 V 

n — 1 j — 1 

1/2 A L''~ 7 J 



1/2 



^Myn) E( n+1 ) d " 3 E 3 

\ p 7 n=l j=l 



-3/2 



1/2 



ft" 1 / 2 . 



A 



1 



(4.21) 



^E(^-«+i)(n+D d - 3 E (*_ B+1)j 

n=l 3=1 v /J 

< ck/i 1 " 7 < c^ 1/2 ~ 7/2 , 



t (kj) x 3 
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and 
(4.22) 



N 



lh-<\ 

h d-iJ2{N-n + l)(n + l) d - 3 



6 



<chlog(l/h), 



^ (N-n+l)j 

n=\ j=l v ,J 

Combining ( [4~18| - ( |4~22| ) with §4~17\ gives |412| . 

We take care of the case of a non-trivial number of zeroes of G" on [0, R] exactly as in 
the d = 2 case. 

□ 



We can now prove Theorem |1.1| 

Proof of Theorem \l.l\ Theorem |1.1| is an direct consequence of Proposition |3TT| Lemma |4.1[ 
and Proposition |4.2| □ 



5. Examples of potentials that satisfy assumption 11.81 

We begin by deriving an integral expression for the scattering angle, E(a). Note that, if 
the potential is non-trapping at energy E, then along a bicharacteristic, the functions p and 
p do not any have simulateous zeros. For if there were such a time, and the value of r at 
this time were ro, then r(t) = tq would be a bicharacteristic, contradicting the non-trapping 
assumption. Hence the zeros of the function E — a 2 /r 2 — V(r) are simple on the region of 



interaction (1.4). Given a fixed bicharacteristic, let r m be the minimum value of r; note 
that r m is a function only of the angular momentum a. We denote the derivative of r m with 
respect to a by r' m (a). By symmetry, this is the unique point at which p = and p > 0, so 
we can divide the bicharacteristic 'in half,' and consider only times when r > r m and p > 0. 
We write 

dip a/r 2 a 
dr p r 2 y/E - a 2 /r 2 - V(r) ' 
By the simplicity of the zeros in the denominator, we can integrate to obtain, for a > 0, 



(5.1) S(a)=7r-2/ dr. 

Jr m r 2 ^E -a 2 /r 2 -V(r) 

(See Figure[l] The change in 9 is the angle between the incoming and outgoing trajectories.) 

Proposition 5.1. Suppose that on the region of interaction 1Z the potential V satisfies, V' < 
(so necessarily V > since V is compactly supported) and r(V' ) 2 + (E — V)(V + rV") > 
with equality only when V — 0. Then 9 Q E(a) < for for a € [0,-R). 

The same is true ifV = cW , where c is sufficiently large and where W{r) — for r > R, 
W(r) > for < r < R and W"(r) is positive and monotone decreasing in some nonempty 
interval [R — e, R) . 



In (5.1 ), set s = r/r„ 



S(a) = 7r — 2 / . d(sr rn ) 

' (sr m ) 2 ^E-a 2 /{{sr m ) 2 -V{sr m ) 

rv 

ds. 



s 2 y/Er 2 n - a 2 js 2 - r 2 % V{sr m ) 
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Differentiating under the integral sign gives 

1 r°° ( 1 



2 Ji \s 2 ^/Er 2 n - a 2 /s 2 - r^V{srm) 

a 2r m r' m - 2a/ s 2 - 2r m V(sr m )r' m - r 2 m sV {sr m )r' n 
2 s 2 (Er 2 n — a 2 / s 2 — r^ n V(sr m )) 3 ^ 2 

(r 2 m {E - V(sr m )) ~ ar m r' m (E - V(sr m ) - \r m sV (sr m )) 



ds 



s 2 (Er^ - a 2 /s 2 - r 2 n V{sr m )f' /2 
Differentiating 

(5.2) E - a 2 /r 2 m - V(r m ) = 

shows ar m r' m = E ^ v{rm) fi rmV ^ rm)y Plugging this in gives 

r(ji,r, t,/„ ^ 2 E-V{sr m )-\r m sV'{sr m ) 



ds 



S'(a)= / [r 2 m (E-V(sr m ))-a 2 



2 KJ Ji V™ ' E-V{r m )-\r m V'{r m ) 

ds 



s 2 (Er 2 ^ - a 2 /s 2 - r^F(sr m )) 3/2 



and using (5.2) again shows that (l/2)<9 Q E(a) is equal to 
E - V(sr m ) - \r m sV'(sr m ) 
s 2 {Er 2 m ~ a 2 Is 2 - r^F(sr m )) 3/2 

E - V(sr m ) E - V(r rl 



r^ds 



E - V(sr m ) - \r m sV'(sr m ) E - V(r m ) - \r m V'{r m ) 
E-V(r) - \rV'{r) 
r m r 2 (E-a 2 /r 2 -V(r)) 3/2 

E-V{r) E - V{r m ) \ 

E-V(r)-\rV'{r) E - V{r rn ) - \r m V'{r m ) ) 

Differentiating the expression (E— V(r))/(E — V(r) — \rV'(r)) with respect to r and using 
V < 0, we see that the integrand is positive if for r m < r < R if 

(5.3) r{V') 2 + {E-V){V +rV") > 0. 

This is exactly the condition in the first paragraph of the proposition. 

To check that the condition in the second paragraph is also sufficient, observe the follow- 
ing. If V(r) > 0, V'(r) < and V"(r) > on some open interval (R — e, R), it follows that 
for r sufficiently close to R, that V + rV" > 0. By picking large enough c > 0, (5.3) will 
hold on the region of interaction 1Z. 



6. Scattering by the disk 

In this section we will prove Theorem 1 1 . 5| from the introduction. We restrict our attention 
to the disk of radius 1, since the phase shifts for the disk of radius R can be obtained from 
those for R = I by a scalling argument. 

Here we use on an ODE analogous to that in ( |I.3| to give a formula for the eigenvalues. 
In fact, for any smooth solution /; to Agd-i/j = l(l + d—2)fi, a straightforward computation 
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shows that 
(6.1) 



„ H l + (d-2)/2( k ) , 
bk(Jl) - ~T^2) 77ZJU 



Z+(d-2)/2 

(i) m 

where the He, are Hankel functions of order v [I] . 
For z > dchnc £ = ((z) by 



(k) 



(6.2) 



3^ 

f(-0 



3/2 



1 VT^ 



dt = log 



l+vr 



z\ 



3/2 



t 



-dt = \f 



z 2 — 1 — cos 



The two formulas for £ are equivalent; the first is more convenient for z < 1 and the second 
for z > 1. 

Lemma 6.1. ^4s fc, / — > oo 

(6.3) 

+ 0((!/- 1 C 3/2 )°°) + 0(^°°), k/v<l, 
+ arg (exp (-|K 3/2 )) + OC^-^" 3 / 2 ) + 0{v-\V 2 ) 

+0(i/-°°), fc/i/ > 1 



arg iff 1 ) (fc) = < 
where v = I + (d— 2)/2, and 



C = C(*/f ) 



(see (|6.2D;. 

In particular, for any p £ (1/3,1), «/ we restrict to the regions k p < v < k — k p or 
v > k + /c p , tten £/ie error terms in (6.3) becomes 



(6.4) 



arg 



— 7T 

2 

-5 
12 



0(kr°°), v>k + kP, 

f + arg (exp (~f K 3/2 )) + 0(fc (1_3p)/2 ), k p < v < k - k p 



We will prove this in a moment. Using Lemma |6.1| and the arguments above, we can now 
prove Theorem 1 1.5| 



Proof of Theorem \1.5\ First we will bound D(£ 1 / k (e)), where £x/fe( e ) is as m (4.2) G(a) = 



\a(,(l/a) 3 / 2 and a — vjk. We apply (4.12), using 



(6.5) 



k < max G < c < oo 

0<Q<1 

p > min G"(a) > c> 0. 

0<a<l 



This means that in (4.12), p(e, 1/k) > c > and k(e, X/k) < c < oo for all e, and thus for 
any < 7 < 1, 



(6.6) 



D(£ 1/k (e))<c(k-i + k-V 2+ ^ 2 ) 



Now we bound the true driscrepancy £\/k by following the proof of Proposition 



4.2 



with 



k = l/h, sp ecific ally using (|47| and C(e)h = 0(k^- 3p ^ 2 ), e = k^ 1 and taking 7 = 1/3 in 
Q gives ( fTl2| . 

□ 
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Proof of Lemma 6.1' We will use the following uniform asymptotic (9.3.37 from |1J. For 
C = C(z) as in mW 



(6.7) 



Hl(vz) = 2er^ 3 



4C \ 1/4 ( Ai (e 2 "/ 3 ^ 2 / 3 C) e 2 "/ 3 Ai' (e 2 ™/ 3 */ 2 / 3 C) 



1-z 2 



I/V3 



^5/3 



6o(0 



x (1 + O fl (i/- a ) + 0(1/-°°)) , 



where z/ —> oo, z £ (0, oo), bo is bounded and positive, and Or denotes a real- valued 
remainder. 

The Airy function Ai(x) satisfies 



Ai(x) ~ i 7 r- 1 /2 :E -i/4 e -2 a = 3/2 /3^ Cfe;r - 



(6.8) 



* -l/2„.l 



fe=0 
oo 



Ai'(x) - — ^-7r~ ' x ' e 



/2„l/4 -2x J ' 2 /3 



5> 

fc=0 



3fc/2 



-3fc/2 



as |x| — > oo with |argx| < ir (chapter 10 of pQ). The Ck and are real, and cq = do = 1. 
We will plug (|6]§ into (JHTTj) with x = e 27r4 / 3 i/ 2 / 3 C- 

(6.9) 

7T /Ai(e 2 "/ 3 ^/ 3 C) e 27ri / 3 Ai'(e 2 "/ 3 i/ 2 / 3 C) \ 
axgif>z) = -- + arg 1 ^ + ^ ^feo(C) + 0(^°°) 



= -- + argAi (x) 



+ (,-~) + .r g (l-(5^fl)) + 0(,-~). 



Consider first the case z < 1. Then £ > 0. Using that arg Ai(x) = — tt/6 + 0(x °°) and 
that Im(e 2 " /3 Ai'(x)/ Ai(x)) = 0(x _0 °), we conclude that 

(6.10) arg #1(1/2;) = -tt/2 + 0(i/ -00 ) + 0(s -00 ). 



Now consider the case z > 1. In (6.9) we take an additional step 



arg If* (i/z) = - - + arg Ai (x) 



°(^) + ^( i -( £ ^S 1 )) +o "' 



7T 

3 



(x^e-^'/^l + OCx- 3 / 2 ) 



(6.11) 



+ 0(i/-°°) + axg 1- 



3 2 ™/ 3 Ai'(x) 



Ai(x)*/ 4 / 3 



arg I x ' e 



1/4 -2f( 3/2 /3 



3 ) +0(x~ 3 / 2 ) 



+ 0(j/- oo ) + 0(x 1 / 2 */- 4 / 3 ). 



This completes the proof of (6.3). 



□ 
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